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I We present a full quantum theory for the dissipative dynamics of an optical cavity in the ultra- 

. strong light-matter coupling regime, in which the vacuum Rabi frequency is a significant fraction 
of the active electronic transition frequency and the anti-resonant terms of the light-matter cou- 

f*H I pling play an important role. In particular, our model can be applied to the case of intersubband 

, transitions in doped semiconductor quantum wells embedded in a microcavity. The coupling of the 

i—i • intracavity photonic mode and of the electronic polarization to the external, frequency- dependent, 

I ' dissipation baths is taken into account by means of quantum Langevin equations in the input-output 

, formalism. In the case of a time-independent vacuum Rabi frequency, exact analytical expressions 
for the operators are obtained, which allows us to characterize the quantum dissipative response 

I I ' of the cavity to an arbitrary initial condition (vacuum, coherent field, thermal excitation). For a 

^ , vacuum input in both the photonic and electronic polarization modes, the ground state of the cavity 

(— I ■ system is a two-mode squeezed vacuum state with a finite population in both photonic and electronic 

' modes. These excitations are however virtual and can not escape from the cavity: for a vacuum 

O I input, a vacuum output is found, without any trace of the intracavity squeezing. For a coherent 

. photonic input the linear optical response spectra (refiectivity, absorption, transmission) have been 
studied, and signatures of the ultra-strong coupling have been identified in the asymmetric and 

M I peculiar anticrossing of the polaritonic eigenmodes. Finally, we have calculated the electrolumines- 

I , cence spectra in the case of an incoherent electronic input: the emission intensity in the ultra-strong 

' O ' coupling regime results significantly enhanced as compared to the case of an isolated quantum well 

C ' without a surrounding cavity. 

In recent years, the study of cavity quantum electrodynamics^^ has been a subject of an intense and fruitful research. 
' In particular, a considerable deal of activity has been devoted to the so-called strong light-matter coupling regime. This 
^ , regime is achieved when the vacuum Rabi frequency of an electronic excitation (i.e., the Rabi frequency associated 
' to electric field vacuum fluctuations) exceeds the frequency broadening of the photonic and electronic excitations. In 
the case of atoms in high-finesse cavities (either in the optical^ or in the microwave^ domains), the strong coupling 
occurs thanks to the strong spatial localization of the electromagnetic field and the very high quality factors (as high 
25^ as 10*) of the cavity mode and the atomic resonance, but the vacuum Rabi frequency ^Ir still remains a tiny fraction 
i: of the active transition frequency uieg ■ In the case of semiconductor planar microcavities^, the quality factors are by 
■ far less impressive (typically not much larger than Q « 10"^), but the strong coupling can be still comfortably obtained 
' thanks to the strong photon confinement and the large electric dipole moment of the active electronic transitions in 
the quantum well. Both these facts cooperate to give a much larger vacuum Rabi frequency. Using inter-subband 
, transitions in doped semiconductor quantum wella^*^*Si2i^, it is even possible to achieve an unprecedented ultra-strong 
I ■ light-matter coupling regime^, in which the vacuum Rabi frequency flfj becomes comparable to the intersubband 
' O [ electronic transition uji2. In this regime, the standard rotating wave approximation is no longer valid, and anti- 
^ ■ resonant terms of the light-matter interaction start playing a significant role. In particular, it has been shown in Ref|^ 
that the ground state of an isolated cavity is a two-mode squeezed vacuum containing a finite number of virtual and 
correlated photons and intersubband excitations. 

In the present paper, a full quantum description of the system is developed, which now explicitely includes the 
coupling to dissipation baths. In the specific case of intersubband excitations in planar semiconductor microcavities 
''^^ considered here, the dissipation baths are not only responsible for damping rates which are typically as large as 5-10 % 
of the transition frequency^i^, but also provide the way of exciting and observing the cavity dynamics. The photonic 
mode is coupled to external world mostly because of the finite reflectivity of the cavity mirrors, while the intersubband 
transition is coupled to other excitations in the semiconductor material, e.g. acoustic and optical phonons, and free 
carriers in levels other than the ones involved in the considered transition. In particular, the coupling to this electronic 
bath allows one to electrically excite the intersubband transitions. 

The theory here developed is based on the so-called input-output formalismiSiiiiiS, in which the dynamics of the 
electronic polarization and cavity photonic fields is described in terms of quantum Langevin equations for the two 
coupled quantum fields. Differently from previous treatments, we have to take into account here all the anti-resonant 
terms of the vacuum Rabi coupling, which forces us for consistency to keep track in an exact way of the frequency 
dependence of the dissipation baths. 
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The paper is organized as follows. In Sec. we introduce the model Hamiltonian for the considered system and for 
the baths. In Sec. ^ the quantum Langevin equations for the intracavity photonic and electronic polarization fields 
are derived. Their exact solution for the intracavity operators as well as for the output operators is given in Sec lIIII 
The case of a vacuum input both for the photonic and polarization fields is analyzed in Sec. IIVI where the properties 
of the ground state of the system arc characterized and the presence of squeezing in the intra-cavity fields is pointed 
out. The question of the observability or not of this squeezing in the output field is addressed in Sec. where we 
predict that no squeezing can be observed in the output modes unless the input is itself squeezed. 

After the brief discussion of single- and double-sided cavities of Sec lVI Al the main features of the linear optical 
spectra (reflectivity, absorption, transmission) under optical excitation are considered in Sec lVI Bl The electrolumi- 
nescence emission spectra under an electronic excitation are studied in Sec lVI 01 where a remarkable enhancement of 
the emission intensity is found as compared to the case of an isolated quantum well in free space. Conclusions and 
perspectives are finally drawn in Sec. IVIII 



I. MODEL HAMILTONIAN 



As schematically shown in FigQ] the Hamiltonian of the present system contains three main blocks: 

H = H,ys + HP,!:,t + H§i,^, (I) 

where Hsys is the Hamiltonian describing the closed cavity system, while Hf^^^f'^ and H^jj^^f^ take into account the 
coupling to respectively the photonic and electronic reservoirs. 

The Hamiltonian Hgys = Hq + Hres + Hanti for the closed cavity system has the typical Hopfield^^ form already 
discussed in detail in Ref. ^ and contains three parts respectively describing the energy of the bare cavity photons 
and of the electronic excitations (Hq), the resonant part of the light-matter interaction (Hres), and the anti-resonant 
terms^i usually neglected in the so-called rotating- wave approximation^ - (Hanti)- 



Ho = ^?iwcav,fc (alak + ^ j + ^ ?iwi2 b^bk , (2) 

k ^ ' k 

k 

Hanti = TT.'^{i^R,k{akb-k- a^k^^_k) + Dk{aka-k + a}ka}_^.)} . (4) 



Translational invariance along the cavity plane has been here assumed. is the creation operator of a cavity photon 

with in-plane wave-vector k and energy ?iti^cav.k- b 
mode of wavevector k of the doped quantum well^: 



with in-plane wave-vector k and energy ?iti^cav.k- &k creation operator of the bright intersubband excitation 



Here, Nqy/ is the number of quantum wells present in the cavity (which are assumed for simplicity to be situated at 

the antinodes of the cavity mode, and therefore identically coupled to the photonic mode), aei = Nei/S is electron 

(j) 

density per unit area in each quantum well and S is the quantization area. The fermionic operator c^^ annihilates an 

electron of in-plane wavevector q from the lowest subband of the j-th quantum well, while C2 q creates an electron 
with wavevector q in the second, excited, subband of the j-th quantum well, kp is the Fermi wavevector of the 
two-dimensional electron gas in each well, and the electronic ground state is written as: 

Nqiv 

1^) = n n ^I'q^ (6) 

j=l |q|<fci? 

in terms of the empty conduction band state \Ocond)- In the following, we shall always restrict ourselves to a weak 
excitation regime: 



^E(^kfck)«^ei, (7) 
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in which the operators 6k are approximately bosonic: 

[bk,bl] =SkM'- (8) 

The NqwNei — 1 states orthogonal to b^\F) are dark excitations (uncoupled to the cavity photon) which do not play 
any role in the optical processes taking place in the system. The first and second electronic subbands have been here 
assumed to be perfectly parallel22,, with a momentum-independent transition frequency lji2. 

^R,k is the vacuum Rabi frequency which quantifies the strength of the light-matter dipole coupling: in the case of 
intersubband transitions in doped quantum wells^, it can become a significant fraction of the intcrsubband transition 
uji2. The explicit expression is 



^R,k = ( i-g- cFei Nqw fi2 sin^ ^ ) : (9) 

Veoo"loi^cav / 

where i^av is the effective length of the cavity mode and 9 is the intracavity photon propagation angle such that 
sin 6* = ck /{uj 12 y/e^). 

Dk is the related term originating from the squared electromagnetic vector potential term in the minimal coupling 
light-matter Hamiltonian^: for a typical quantum well potential, one has D)^ ~ j,/tt'i2. 

The environment of the open cavity system is modeled by two baths of excitations, associated to respectively the 
photonic and the electronic degrees of freedom. In this paper, we shall focus our attention on cavity configurations 
in which the photonic mode is coupled to the external electromagnetic field via the finite transmittivity of the planar 
mirrors enclosing the cavity, as recently don o'^i^i^ . In a planar geometry, the coupling of the cavity photon to the 
extra-cavity electromagnetic modes is well described by the Hamiltonian: 

Hbath ^ I dqY^ ^^IX ( + ^j+ih 1 dqJ2 {^t^ "g.k 4 - K^lL* "U'^k) , (10) 

k ^ ^ k 

where is the frequency of an extra-cavity photon with in-plane wavevector k and wavevector q in the orthogonal 

direction and is the corresponding creation operator, obeying the commutation rule [a^.k, aj/ k'] = 9')'^k.k'- 
The coupling between the cavity and the extracavity radiation fields is quantified by the tunneling matrix element 
'^q'k through the cavity mirror, whose value depends on the specific mirror structure and can be calculated by solving 
the classical Maxwell equations. 

The Hamiltonian H1U|I takes into account a single radiative bath, coupled to the cavity through one of its mirrors. 
For the sake of simplicity, we shall start our analysis from this simplified case, which is then extended in Sec lVI Al 
to the more realistic case of two photonic baths coupled to the cavity through the front and the back mirrors. All 
this discussion is performed by neglecting the anti-resonant terms in Hfl^°^ whose effect is quantitatively small; the 
general case including would not however pose any additional problem, and is considered and solved in Appendix 1X1 

Concerning the bath coupled to the electronic transition, the modelling is not as straightforward as in the photonic 
case. In the case of intersubband transitions in doped semiconductor quantum wellsi^ii^, the damping and dccohercnce 
of the electronic transition at the frequency uji2 is due to the interplay of different processes, such as the interaction 
with crystal phonons (optical and acoustical), the scattering with impurities and with free carriers in levels other 
from the ones involved in the considered electronic transitiorki&. As our purpose is to study the role of the electronic 
losses and decoherence on the peculiar quantum optical properties of cavities in the ultra-strong coupling regime, we 
do not attempt here to address the specific nature of the decoherence channels, but we will rather model them as a 
phenomenological bath of harmonic excitations in the spirit of the Caldeira-Leggett modelii. As for the photonic bath, 
the harmonic oscillators modes coupled to the intersubband excitation of wavevector k are labelled by a continuous 
index q and have frequency w^'^.. They affect the electronic intersubband transition according to the Hamiltonian: 

Here, the bath operators /3q,k satisfy the harmonic oscillator commutation rule [/3g,k,/3^/ k'] = ~ 9')'^k,k'- k^'^. 
are the matrix elements quantifying the coupling to the electronic polarization. Inclusion of the anti-resonant terms 
in H^l^f^ is straightforward, as discussed in Appendix 1X1 In practice, ti;|'^ and At^'j^ can be chosen at will so as to 
quantitatively reproduce the properties of a specific system, in particular the frequency-dependent damping rate of 
the electronic excitations. A recapitulative scheme of the investigated model in depicted in Fig. ^ For a more detailed 
description of the specific system of semiconductor microcavities embedding intersubband transitions, the reader can 
refer to Refs. and0. 
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FIG. 1: A sketch of the investigated model with the corresponding operators and frequencies. The system is represented 
by a planar cavity mode embedding a doped multiple quantum well structure. The cavity photon field is coupled to the 
quantum well electronic polarization associated to a resonant intersubband transition in the doped quantum wells. The 
Hopfield Hamiltonian for the cavity system includes the anti-resonant light-matter interaction terms, which are significant in 
the ultra-strong coupling regime^ (vacuum Rabi frequency Q.R^k comparable to the transition frequency lo\2)- The cavity mode 
is coupled to the extracavity electromagnetic field. The electronic polarization is coupled to a bath of electronic excitations. 



II. QUANTUM DYNAMICAL EQUATIONS 
A. Quantum Langevin equations 

The equation of motion for the extra-cavity photon operator in Heisenberg representation reads: 



dt h 

and its solution can be formally written as 



ph * 
k 



at e 9'''^ ' 



ak(i') 



(11) 



(12) 



to being the initial time. Inserting these formulas into the evolution equation for the cavity photon amplitude, one 
finds: 



~dt 



^9 KX "S'k 



= - J[ak, H,y,] +Jdq e-^'-^*"*") - j dq |fc,,k|' dt' e-K,t(*-*') a^{t') . (13) 



Using the standard definition 



*9,k 



"g,k(io)e 



for the input fields at ~* — cxd, one can cast in the form of a quantum Langevin equation: 

= ^^y"^ ^ I rcav,k(i - t') ak{t') + Fcav,k(i) , 



(14) 



(15) 
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FIG. 2: A recapitulative sketch of the input-output scheme. The initial condition on the extracavity photon operators and the 
electronic bath operators at t = — oo represents the input. The intracavity dynamics for the cavity mode and electronic polar- 
ization field is described by quantum Langevin equations including the antiresonant terms of the light-matter interaction. The 
radiative and non-radiative baths produce a complex frequency-dependent damping as well as frequency-dependent Langevin 
forces. The bath operators at f = -\-oo represent the radiative and non-radiative output. 



where the (causal) damping memory kernel is given by 



(r)=e(r) / dq\nfj'e~<^^ 



and the fluctuating Langevin force is represented by the operator 



Fcav,k(t) ^ j dq ^^\e-'<"-'o^^^ 



Analogously, one obtains a quantum Langevin equation for the electronic polarization field 

96k i 



^ [bk, Hsys] - / dt'Ti2,),{t - t') fek(i') + Fi2.k(i) 



in terms of the memory kernel 



ri2,k(r) = e(r) / dq Kkl'e^*'^--^ , 



and the Langevin force 

Fi2M{t) ^ j dq 4^e-^<-'(3\\. 
The input operators /3*\ for the electronic bath are here defined in the same way as the photonic ones H14|l. 



(16) 



(17) 



(18) 



(19) 



(20) 



B. Input-output relations 

The extra-cavity asymptotic output operators &t t = +oo can be related to the input operators aX to — —oo and the 
cavity photon ones through a linear relationship. Taking — > — oo and t — s- -foo in Eq. p2|l . we obtain the formula 



„ out „ in „P'»* ~ ,ph \ 



(21) 



where a\^{u)) is the Fourier transform of ak(^)- An analogous expression holds for the electronic bath operators: 



pout Qin .^eL* u ( , .et \ 
Pg,k = Pq.k - ^q,k Ok(Wg,k) 



q,k; 



(22) 
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C. Equations in frequency space 

In the present case of a time- independent Hsys , the quantum Langcvin equations (|15|l and H18() are most conveniently 
solved in the frequency space. By Fourier transforming them and their hermitian conjugatesiS, we get the following 
equation 



Ml 













+ i 




V 



(23) 



where the 4x4 matrix 



V 



/ t^cav.fc + '2Dk - LU - ircav,k(t^) i^R,k 2L»fc -i^R,k 

-i^R,k uj 12 - uj - iT 12, -i^R,k 

-2L)fe ~i^R,k -Wcav,fc - 2L)fc -UJ - if cav.k("^) i^R,k 

-iflR,k -i^R,k -Wi2 - w - if|2 k(^^) / 

(24) 

is the Hopfield^*^ matrix of our system, additioned by the terms accounting for the damping and for the frequency 
shift produced by the coupling to the baths. The complex self-energy shift of the cavity photon due to the coupling 
to the photonic bath is indeed 



rcav,k(w) ^ I dq TT ~ ^g,k) + t'P j dq 



K 



g,kl 



U! - UJnM 



(25) 



Its real part represents the frequency-dependent radiative damping of the cavity mode due to radiative losses. This, 
in agreement with the usual Fermi golden rule 

5R[fcav,k(c.)]=^|fcfklVf(^), (26) 

where p^{uj) = [d^q'^q dq] ^ is the photonic density of states at in-plane wavevector k and q is the resonant 

wavevector such that a;?'j^ = oj. The imaginary part 3[fcav,k(w)] accounts instead for the corresponding Lamb shifli^: 

for a given 3?[rcav.k(^)], the Lamb shift 3[rcav,k('j-')] is in fact univocally fixed by the Kramers-Kronig causality 
relationships 



5{fcav,k(t^)} = —V 



duj 



^3fi{fcav,k(c<^')} 



An analogous formula holds for the electronic excitation counterpart: 

fi2,k(^) = fdq |4^kl'^'5(a;-<k)+«^ / dq 



^el 12 
^g,kl 

el 

- ^q,k 



(27) 



(28) 



The real part of ri2.k('^) represents the frequency-dependent broadening of the electronic transition, of non-radiative 
origin. As all real excitations of the considered baths have by definition positive frequency i^^k''' > 0, one has: 

3?{f eav,k(c^ < 0)} = 3fi{f i2,k(w < 0)} = . (29) 

In frequency space, the Langevin forces have the form 

Fcav,k(u;) ^ [dq 2n5ic. - co^i) o^^\ = 2n^fpi'iu;) a^^ (30) 



Fi2Mi^) - / dq' 27t6{cu - ojp^^) (3^^ = 2^^'^;, p^^{lo) (i^^^, 



el el I 



(31) 



where q and g' in the right-hand sides are such that i^^k ~ k ~ analogy with H29fl . one therefore has 

Fcav.k(w < 0) = Fi2,k(w < 0) = . (32) 

A recapitulative sketch of the input-output framework here developed is drawn in Fig. [21 
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III. EXACT SOLUTIONS FOR THE OPERATORS 



The equation H23|) for the Fourier space cavity field operators is immediately solved by matrix inversion; 



/ ^l;av,k(^) \ 

V ^iV-k(-^) ; 



(33) 



where 



For ui > 0, this can be simplified as 



ak(w > 0) = t/ii(k,a;) Fcav.k(w) +t/i2(k,a;) Fi2,k(w) , 
6k(w > 0) ^ ^21 (k, a;) Fcav,k(w) +^22(k, a;) Fi2,k(w)- 



(34) 

(35) 
(36) 



Using the expressions Eqs. (|21I22(I for the output field and the expressions H3()I31|I for the Langevin forces, one finally 
gets to the input-output relation: 



^ _ f njiiKu^) z^i2(k,c^) 



fJOUt 

Mg',k 



Z^2l(k,t^) Z^22(k,t^) 



a 



in 

V,k 



(37) 



This is a linear relation linking the output operators of the photonic and electronic bath modes to the corresponding 
input ones; as expected by the time-invariance of the Hamiltonian under examination, the pair of photonic and 
electronic modes involved in H37I) share by construction the same frequency 



ph 



uj„, t . The matrix elements of U 

q ,K 



are 



L(ii{k,Lu) = 1 - 23fi{fcav,k(w)}aii(k,w) , 

ph * 

Z^i2(k,c^) = -25R{fi2,k(a;)}^ai2(k,c^) , 

Z:?2l(k,C^) --23fi{feav,k(w)}^^2l(k,C^) , 

%,k 

Z^22(k,C^) = l-23fi{fi2^k(c^)}^22(k,C^) . 



(38) 
(39) 

(40) 
(41) 



The unitarity of the matrix U(\i, uj) (that has been numerically checked) guarantees that the usual Bose commutation 
rules hold for the output operators a°"^ and /?° k*' namely [a^^, ^^"k'] = Hi^q') '^k,k' and [0°'^^ ^ Pq^U] = S{q—q') 5k, k'- 
^From a more physical point of view, the unitarity of U (k, lo) means that the total energy is conserved during a 
scattering process when energy is sent onto the cavity as radiation or as electronic energy, and both the emerging 
light and the electronic absorption are taken into account: 



out t otit I QOut\ QOut\ / „in i t om \ 

"9,k "g.k + Pg',k Pq'M ) = \",,k "<;,k + Pg'.k Pq' M ) 



(42) 



It is interesting and an important check of consistency to verify that the usual commutation rules [ak(t), aj^., (t)] = 
<Jk,k' hold for the cavity (as well as for the electronic polarization) operators. This has been verified on our specific 
model in the following way. By definition, one has: 



Inserting here and noticing that 



-iiut Auj t r~ 



e-'*[ak(a;),4,(c.')] 



(43) 



[Fcav,k(w),Felv.k'(^')] - 2^-5(60 - iu')2!fi[t. 



k.k' 



(44) 
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FIG. 3: Quantum properties of the intra-cavity ground state for input fields in tlie vacuum state. In this case, the output is 
also in the vacuum state, which implies that the cavity photons present in the cavity are purely virtual ones and can not escape 
the cavity. Left panel: cavity photon number (aj^ak) (solid line) and the real part of the anomalous expectation value (akfl-k) 
(dashed line) as a function of the normalized vacuum Rabi frequency flii/uji2. Right panel: corresponding cavity photon 
correlation Ck,-k (solid line) and the normalized difference noise Fdiff (dashed line). oJcav.fe = (^12, and the damping rates are 
taken as frequency-independent ones. The calculations have been performed with 5R{rcav,k(<jJ > 0)} = fcav,fc ~ 3^{rcav,k('^ > 
0)} = Hav.fc = 0.04t^i2 and K{ri2,fc(a; > 0)} = fi2,fc = O.O81.J12. 



one can finally write 



2^ 



|gi3(k,c^)|2 5R[fcav,k(-^)]- I^i4(k,u;)|'3?[fi2,k(-^)] , (45) 



which has been numerically checked to give 1 as expected. This fact critically depends on the consistent inclusion of 
the real and imaginary parts of rcav,k(w) and ri2,k(w) satisfying the causality relation (|27() . 

Note finally that the presence of the anti-resonant terms in Hsys implies that the cavity operators a]s_{uj) and &k('-^) 
have non- vanishing values also for negative frequencies: 



ak(ij-' < 0) 



^i3(k,^) Fl^_^i-cj) + gii{k,u) 



6k(c^ < 0) = fe(k,w) Fl^j-oj) + g24iK^) ^iV-k(-^) • 



(46) 
(47) 



IV. THE CASE OF A VACUUM INPUT: PROPERTIES OF THE INTRACAVITY QUANTUM 

GROUND STATE 



Once we have obtained explicit expressions for the output operators in terms of the input ones, we can apply this 
formalism to study the quantum dissipative response of the cavity system for different initial states. 
We shall start from the simplest case of a bath initially in its vacuum state |0) such that: 

<k|0) = /3;?k|0) = (48) 

for every k and q. As expected on physical grounds, no excitations are created, and the bath remains in its vacuum 
state. It is in fact a straightforward consequence of H37() that 

<k|0)=/3rklO)-0- (49) 

On the other hand, the cavity system itself is not in the usual vacuum state, and a finite number of photons is present 
in this anomalous ground state. This can be calculated by means of H33|l . which relates the in-cavity field at to the 
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input bath operators -Fcav,i2 and F^^^ ^2'- 

= -/ da;|^i3(k,-u;)p5R{rcav,kH} + |ai4(k,-c^)p5R{ri2.kH}. (50) 
I" Jo 

As the input bath is here taken in the vacuum state Fcav,i2 |0) = 0, only the terms involving F}.^^ -^2 contribute to the 
result H50|) . The presence of the anti-resonant light-matter coupling Hanti is here crucial, as it is responsible for the 
non-vanishing value of the matrix elements ^13 and t/14 connecting creation and destruction operators. 

Remarkably, the ground state shows finite anomalous correlations between modes with opposite wavevectors: 

(ak(t)a_k(0> = TTT^ / duo / dio' (4(c.)ak(c.')) e'^^+^ = 
[^^) J-00 J-00 ^ ' 

= - f dw^n(k,w)^i3(-k,-w)3?{feav,kM}+fe(k,C^)^14(-k,-C^)3?{fi2,kH} , (51) 
Jo 

which suggest the ground state to be a sort of squeezed vacuum. 

The predictions (|50|l and (|51|l for respectively the number of photons (solid line) and the anomalous correlations 
(dashed line) have been plotted in the left panel of Fig. |2| as a function of the normalized vacuum Rabi frequency 
flji/LUi2 for the resonant case Wcav,fc = 1^12- For simplicity, in Figl^as well as in all the following figures, a frequency- 
independent damping of the cavity and electronic excitation modes has been used, i.e., 3fi{rcav,k(^^ > 0)} = Fcav,*; 
and 5R{ri2,k(w > 0)} = f 12,^. Note, however, that our theory is able to describe the system under consideration for 
colored baths with an arbitrary frequency dependence of rcav,k(w) and ri2,k('^). 

It is interesting to check that in the limit Fcav.k, ri2,k ^ the values of (|50() and H51|l are in quantitative agreement 
with the result of a direct diagonalization of the isolated cavity Hamiltonian Hgys as done in Ref 9. Up to moderate 
values of the broadening values, the difference from the isolated cavity preduction remains small, in particular the 
total number of virtual photons is only slightly changed. 

Insight on the structure of the ground state can be obtained by looking at other observables which are most sensitive 
to two-mode squeezing effects. In the right panel of Fig. |2| we have shown the results for the photonic correlation 
Ck,-k (solid line) between the noise amplitudes of the k and — k-modes. As usual, the amplitude quadrature operator 
of the k-mode is defined as = Ok + aj^. In the ground state, this correlation is 

_^ ^ ((^k-(lk))(^-k-(l-k))) ^ 2jR{(aka-k)} ^ ^^32) 
'((^k-(^k))2)((X_k-(^-k))2) l + 2(4«k) 

where we have used the fact that in the ground state (^k) — (-^-k) = 0. As shown in right panel of Fig. O the 
correlation tends to when Qji/uji2 0, i.e., the normal vacuum is reobtained in the weak-coupling limit. 

In the right panel of Fig. |21 the normalized variance of the difference between the amplitude quadratures of two 
correlated modes is shown, namely 

Fdiff = i ((Ik - ^-k)') = 1 + 2(a[ak) - 23fi{(aka-k)} • (53) 

Fluctuations below the shot noise level correspond to F^iff < 1. As shown in the figure, F^is tends to 1 for vanishing 
flji/L0i2. When entering the ultra-strong coupling regime, -fdiff decreases, becoming smaller than 1 and the two-mode 
squeezing of the quantum ground state becomes more and more significant. 

It is important to stress that the discussion in the present section involves the state of the intra-cavity field, and 
that the virtual photonic and electronic excitations present in the cavity ground state are trapped inside the cavity. 
For an input field in the vacuum state, the output is in fact always in the vacuum state as well and no radiation will 
be emitted, in agreement with energy consevation requirements. On the other hand, if the properties of the cavity 
system were modulated in time in a non-adiabatic way, a creation of real excitations is possible inside the cavity, 
which in turn leads to a finite emission of quantum vacuum radiation^ analogous to what is predicted to happen in the 
so-called dynamical Casimir cSeciS^^. A complete study of these issues, and in particular of the emission intensity, 
is presently the subject of detailed investigations^Si. 



10 



V. EXTRACAVITY TWO-MODE SQUEEZING FROM A NON-SQUEEZED OPTICAL INPUT ? 

An interesting question which has been raised in the hterature of interband excitonic transitions^LSSiSS is whether 
any quantum optical squeezingi^ can be observed as a consequence of the finite anomalous correlation shown by 
the polariton vacuumi^. In our specific case, the question we ask is whether the finite anomalous correlations H51|) 
shown by the intracavity photon field can be observed as a squeezing of the output field emitted by the cavity. This, 
obviously, in the absence of any squeezing in the input field. 

The most general non-squeezed optical input state is identified by the condition: 

(«",k')=0, (54) 

which is well satisfied by a thermal incident state, as well as by a coherent incident field^^. The electronic input is 
assumed to be in a thermal state and to have no correlations with the optical input. 

As Eq. (|37|l relates the annihilation operators of the output to the annihilation operators of the input without 
involving the creation ones, one has: 

(alt) = Un{K^=u^',i) («Sc) (55) 
«k c^Tm') - ^ii(k,c. = Lo^^%)Un{K^^ <"k') «k^7M') ' (56) 
from which it is immediate to prove that no squeezing is present in the the output field either: 

(«rk «rk') - «%) = 0- (57) 

This answers our initial question: no squeezing can be ever observed in the output unless the input field is itself 
squeezed, or the properties of the cavity are modulated in time. 



VI. OPTICAL SPECTRA 



So far, the discussion has been limited to the case of a single photonic bath coupled to the cavity. This model is 
sufficient to describe a single-sided cavity in which the cavity mode is coupled to the radiative modes through only one 
of its mirrors, while the other one is supposed to be perfectly reflecting. In order to obtain quantitative predictions for 
quantities of actual experimental interest like reflection, transmission and absorption spectra, one has to generalize 
the model adding a second photonic bath so as to include the radiation emitted by the cavity through its back mirror, 
as sketched in Fig^ 



A. Double-sided cavity: general theory 

The total Hamiltonian must now include a second photonic bath associated to the radiation in the half-space behind 
the cavity, which is coupled to the cavity mode through the back mirror. This can be done by simply including another 
term analogous to l|10|l : 

Hfath ^ J (^lYl ^^qM {^^'Ik^'q.M + \^+^^jdqYl [^t^ «lc - '^^IL * "g!k • (58) 

The new Hopfield matrix J^^°^bie obtained by A^k,t^ in H24() by simply replacing: 

f cav,k('^) f Lv,k('^) = f cav,k(t^) + f cav,k('^) > (59) 

where F^^^^ k(^) complex linewidth of the cavity mode due to the finite transmittivity of the second mirror, 

defined in a way analogous to the definition H25(l of Fcav,k(w). Correspondingly, a new Langevin force -f^cav k('') 
analogous to H17() has to be included, which means to replace in H33|l 

i^cav,k(w) Fl,,kH = Fcav,k(^) + Fc'av,k(^) • (60) 

Following the same steps as in the case of the one-sided cavity, the solutions for the output operators are 

.on* X / ^™ 
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where U'^""'''^'^ {k, is a 3 x 3 unitary matrix. Similarly to the single-sided case of H38I41|I . its matrix elements can be 
written as: 

Ufr'^'i^^) = S,i ~ 25R{fi,k(a;)} 1^^^ gfr'^'iKcu) . (62) 

Here, a shorthand notation has been used: for j — {1,2,3}, the quantities Tj,]s_{uj) respectively mean Fcav.kC'^), 
fi2,k(w), and ^caviii'^)- Analogously, respectively mean k^\^, k^^^^, and k'^'^ and all the three are evaluated at 
wavevector values such that = w^f.k = ^- (k, w) is defined as ^'^°"''''=(k, w) = -i [MdouUe^-^ . 

In the single-sided cavity limit, one has Uf^"''^''(k,uj) = Ujk{Ki^) and Ufp'>^'= (k, cu) = Z^|°"''''='(k, w) = for j,k = 
{1,2}, andUip'''%k,uj) = 1. 



B. Linear optical spectra 

The results of the previous subsection can be used to obtain quantitative predictions for the optical properties of the 
cavity, namely its reflection, absorption and transmission spectra as a function of the frequency uj of the incident light. 
Here, we assume that no input other than optical is present, i.e., = 0, and that the coherent radiation 

incides onto the cavity from the half-space in front of it, i.e. {a^^ ci™k) > 0, while {a^^'' a'^^) = 0. 

In the considered geometry, reflection is described by the output operator a°^, so that the reflectivity is equal to 

7^k(c^) = )% ^' = Kr''^(k,c.)p . (63) 

The transmission through the second mirror is described by the output operator a^"^*- Therefore, the transmittivity 
reads 

^i'^^) = \% = li^'r^^Ku^r ■ (64) 

("g,k "™k) 

Note that 7k (w) = in the single-sided case. The presence of the electronic bath implies that the incident light can 
be absorbed into electronic energy. The corresponding absorption coefHcient of the microcavity system is: 

Au{uj)^\U^r"'iKu;)\^ ■ (65) 
As expected, the total energy is conserved, i.e. 7?.k(w) + ^k('jj) + = 1, as one can verify from the unitarity of 

In Fig. 01 we show an example of reflectivity, transmission, and absorption spectra for different values of the 
normalized detuning S = (wcav.fc — '^12)/'^12, which is varied from —0.5 to 0.5 by steps equal to 0.1. The linear optical 
spectra have resonances corresponding to the cavity polariton eigenmodes. For the large vacuum Rabi frequency here 
considered {rin. = 0.4cji2), it is apparent that the optical spectra show an anticrossing of the polariton eigenmodes, 
which is is strongly asymmetric, as it was anticipated in the case of a closed cavity systemSi without dissipation. This 
can be also seen in Fig. [S] showing the reflection, transmission, and absorption spectra for zero detuning (wcav.fc = 1^12) 
and increasing values of flfi/u}i2 (from to 0.5 with steps of 0.05). Remarkably, note that for an incident frequency 
close to resonance with the polariton eigenmodes, a significative fraction of the incident energy goes is dissipated in 
the electronic bath as absorption. 



C. Electroluminescence 



Another physical quantity which can be successfully studied by the present theory is the intensity of the light 
emitted by the cavity when this is incoherently excited from the electronic bath in a so-called electroluminescence 
experiment. Electrically excited intersubband transitions are in fact playing an important role as light sources in 
the mid and in the far infrared region, but they still suffer from a poor quantum efhciency of radiative emission as 
compared to non-radiative losses. A way of enhancing the emission efficiency would therefore be of great interest. 
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FIG. 4: Reflectivity TZu,{ijj) (left), absorption Au.{uj) (center) and transmission Ti^{uj) (right) spectra as a function of the 
normalized photon frequency uj /u!i2 for different values of the normalized detuning S = (ojcav.fc — (jJi2)/uii2 (from S — —0.5 to 
5 = +0.5 by steps of 0.1). The different curves are offset for clarity. The bottom curve corresponds to 5 = —0.5, while the top 
one has been obtained for 5 = 0.5. Parameters: 57_R.fc/aji2 = 0.4. Broadening parameters as in Fig. |3 




FIG. 5: Reflectivity TZ]^{uj) (left), absorption ^k(i^) (center) and transmittivity Tk{u>) (right) spectra as a function of the 
normalized photon frequency uj/u)i2 for different values of the normalized vacuum Rabi frequency Qii/u!i2 for the resonant 
case ajcav,fc ~ ^12- The different curves are offset for clarity. The bottom curve corresponds to iiu = 0, while the top curve 
corresponds to fin — 0.5a;i2 {Qr/uji2 is increased by steps of 0.05). Same broadening parameters as in Fig0] 



Consider a purely electronic, incoherent input, distributed among the different q modes: 

(Ck^ = I!L, > , (66) 

= (^Z;^^^M)=0- (67) 

The intensity of the spontaneously emitted light through the front mirror (defined as the mean number of photons in 
the q mode) is proportional to the quantity 

C^^ jdq (a™*'t ^out^ ^P/.(^) \atr'"'{k,u;)\'l!Uu;). (68) 
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FIG. 6: A logarithmic scale plot of the dimensionless quantity |Wf2"''''^(k, oj) (proportional to the electroluminescence spec- 
trum) as a function of the normalized frequency Lojiovi in a weak-coupling regime ^rI^vi = 0.001 (dashed-line) and in a 
strong-coupling one Q,r/u\2 = 0.3 (solid line). Other parameters as in Fig. |^ The small increase at very low uj is, & con- 
sequence of the specific choice of constant values for 5R{rcav,k('i' > 0)} and 5R{ri2,k(<^ > 0)}, producing a Kramers-Kronig 
singularity of the imaginary part (Lamb shift) at a; = 0. 




FIG. 7: Solid line: the effective spontaneous emission rate 7^ (in units of 0^12) as a function of the normalized vacuum Rabi 
frequency Q,ii^k/uii2 for ojcav.fe ~ (j->i2- Dashed-line: the Purcell rate 2^^^ ^ri2 ' Horizontal line: the strong coupling limit 
I A/rSravTw. Parameters: fcav,k = fcav.k = 0.04tJi2 and f i2,k = 0.08tJi2. 



The dimensionless quantity \Uf2^'''-'^ (k, a;)^ is represented in Fig. Elfor a weak-coupling situation (dashed-line) and for 
a strong coupling case (solid line), where it exhibits two polaritonic resonances. Suppose that in the spectral region 
where |i7f2"''''^(k, w)p is significant we can roughly neglect the frequency-dependence of the electronic excitation 
(J|^c(a;) « /|ic) and of the extracavity photon density of state {p^{uj) « Pi!^)- With these approximations, the 
luminescence intensity is given by the simplified expression 

-Ck « /|L7f , (69) 

where the effective rate of the luminescence is 

7f = Jdu; \Utr'''{k,uj)\\ (70) 

If we further neglect the frequency-dependence of the broadening, (i.e. 3fJ{rcav,k(^ > 0)} « f ^^^^ ^, ^{Ti2^k{Lo > 
0)} « ri2,k) and consider the resonant case Wcav.k = 1^12, the results can be expressed by simple analytical expressions. 
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In the weak-coupling regime Jljj. j, <C ^^av k^iz-ki it can be accurately approximated by a Purcell-likeaiaai law: 

7k -^fvr — fr —f ' ^'^) 

^ cav,k ^ cav,k ^ 12, k 

where we recall that f ^^^^ = f cav,k + ^cav k the total cavity broadening (due to the front and back mirror). In 
contrast, in the ultra-strong coupling regime, saturates around the value 



l^"" ^ ^V^Lm^i^m (72) 

The numerical dependence of 7^*^ as a function of ^R.k./^i2 is shown in Fig. |7| showing the Purcell-like behavior 
(|71|l in the weak-coupling limit and the saturation value (|72f) in the very strong coupling limit. The slight decrease at 
large f2i?,k is due to the specific shape chosen for the real parts of the damping kernels rcav.k(^) and ri2.k('-^), which 
are constant for a; > and vanishing for lu < 0, with a jump at w = 0. Due to the Kramers-Kronig relationship, the 
imaginary part (the Lamb shift) has a singularity at w = 0. 



D. Comparison with free space case: enhancement of electroluminescence 

It is interesting to compare the predictions for the ultra-strong coupling regime to what occurs for the isolated 
quantum well in the absence of the surrounding microcavity. In order for the comparison of the electroluminescence 
rates to be fair, exactly the same model (|66f) has to be used in both cases for the excitation of the intersubband 
transition by the electronic bath. A huge Hilbert space is in fact available for the electronic excitations of the 
present many-electron system, and the spontaneous emission rate dramatically depends on the specific state under 
consideration In our model, the only bright states are the ones created by the action of the electronic polarization field 
creation operators b^, while the much larger number of other states remain dark. This is in stark contrast with what 
happens in a two-level atom, where the presence of a single excited state makes the spontaneous emission rate to be an 
univocally defined quantity. In particular, an explicit calculation of the quantum efficiency of the electroluminescence 
process with a realistic model of the electronic injection process would require a more refined model taking into 
account the energy that flows into all the dark excited states, orthogonal to the bright mode. For this reason, we 
do not attempt to estimate it in the present paper, but we rather focus our attention on comparing the predictions 
for the emission intensities that are obtained in the two cases using the very same model for the electronic injection 
mechanism. 

In the absence of the surrounding microcavity, no coupling to any cavity mode is present, and the intersubband 
transition is directly coupled to free-space radiative modes. For the frequency-flat excitation of the electronic bath 
the emission intensity can be analytically calculated to be: 

■^k ~ f^k ^exc /bright, k ' V'"^-' 

where 7^]^!^^ ^ is spontaneous emission rate of the bright intersubband excitation of wavevector k when the quantum 
well is embedded in a bulk material of refractive index tQw without any surrounding cavity. 

For an electronic surface density in the quantum well equal to Uei , the free-space spontaneous emission rate of the 
bright excitation mode is calculated by applying the Fermi's golden rule, giving the result 

QW _ 1 NqwCTei „ sin^ 6* 

^''"'5^'''' " 2yeQT^ c 2mo cosO ^ ' 

where e is the electron charge, rriQ is the free electron mass, c the speed of light and Nqw the number of identical 
quantum wells. The propagation angle 9 of the emitted photon is given by sin 9 = ck / {uj 12 y/eqw) ■ The oscillator 
strength /12 of the intersubband transition is written in terms of the electric dipole matrix element Z12 asS^: 

2mo CJ12 

^ ■ V'&j 

Note that the bright state \ F) is the excited state with the largest spontaneous emission rate. In particular, 7^|^ijt k 
increases with the density of the two-dimensional electron gas and does not depend on the emission frequency 0^12. It 
is interesting to compare 7^|^ht k '^ith the spontaneous emission rate of a different excited state of the two-dimensional 
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electron gas, namely C2l^_^_^cfl^\F) . In this case, the intersubband excitation does not coherently involve any longer 
all the electrons in the lower subband, but only the one initially at k, so that the corresponding spontaneous emission 
rate is suppressed by a factor equal to the total number Nqy/N^i of electrons. Note that an atomic-like spontaneous 
emission rate^ (i-e., proportional to 212^^2 and very weak for long wavelength transitions in the far infrared) is 
achieved only for a state like c'2^2\Qcond) (one electron in the second subband, none in the first one). In contrast, the 

strong radiative properties of the bright state b\\F > with spontaneous emission rate 7bHght,k occur in presence of a 
dense electron g considered in this paper. 

For a quantum well with very high barriers /12 does not depend on the quantum well thickness (i.e. on L1J12) and 
is approximately given by /12 — 0.96 mo/m*, m* being the effective electronic mass in the semiconductor-^. Using 
GaAs parameters {tQw — 13.5, /12 — 14), with an electron density Cez = 5 • lO^^cm^^ and Nqw = 10, the radiative 
linewith of the bright state is approximately ?i7^.ight k — 0.04 meV. 

Comparing this result with the one (|72|l for the ultra-strong coupling regime, the cavity-induced enhancement r/k 
of the spontaneous emission of the quantum well is given at resonance — 1^12 by 



- -Ck 7k*^ \/^"IiivM 

^ ~ nQW ~ QW ~ O^rad • ^ ' 

^k 7bright,k 'bright, k 

For an intersubband transition for which 71^^12 ~ 100 meV, typical values for the non-radiative broadening and the 
cavity mode broadening are ?iri2,k, ^Fcav.k ~ lOmeV. The enhancement is then as large as 77k ~ 100. 

The fact that ifT^ holds only in the strong couphng regime wf ^^v k^ i2,k ^ ^^fl,k imposes that: 



^1- « -Sf^ > (77) 

' bright, k 

which means that the enhancement can not made arbitrary large by simply choosing larger linewidths. 

As a final remark, note that the expression (|75|) has been obtained under the resonant condition Wcav.k = W12 and 
is therefore expected to hold in a cone of wavevectors k around the resonant wavevector kres such as Wcav,fcres = '^12- 
On the other hand, when the cavity mode resonance is detuned from the intersubband transition resonance of an 
amount larger than the linewidths, the enhancement with respect to free space case disappears. 



VII. CONCLUSIONS AND PERSPECTIVES 



In conclusion, we have presented a full quantum theory of cavities in the ultra-strong light- matter coupling regime 
(i.e. when the vacuum Rabi frequency is comparable to the active electronic transition frequency), including the 
coupling of the cavity system to dissipative baths. In the case of a time-independent cavity properties, we have 
solved exactly the quantum Langevin equations for the intracavity operators within an input-output approach and 
we have determined analytically the output operators, allowing us to determine the response of the cavity system to 
an arbitrary input. 

In the case of a vacuum input for the photonic and electronic polarization fields, we have characterized the properties 
of the ground state of the cavity system: due to the anti-resonant terms of the light-matter interaction, it turns out 
to be a two-mode squeezed vacuum, whose properties are weakly renormalized by the interaction with the photonic 
and electronic baths. In particular, we have checked that the photons and electronic excitations present in the ground 
state are purely virtual ones, and can not escape from the cavity: if the input is in the vacuum state, the output is 
itself in the vacuum state and no radiation is emitted by the system. Furthermore, we have explicitly shown that if 
no anomalous correlations are present in the input beam, no correlations will be present in the output either, so that 
no extra-cavity squeezing is observable unless the input is itself squeezed. 

On the other hand, the anomalous properties due to the antiresonant terms of the large vacuum Rabi coupling 
show up as a peculiar asymmetric anticrossing of the polariton excitation branches that can be easily observed in the 
optical reflection, transmission, and absorption spectra under coherent light excitation. 

Finally, the input-output formalism has been applied to the study of the electroluminescence emission intensity 
in the case of an electronic excitation: the use of a microcavity surrounding the quantum well provides a significant 
enhancement of the emission performances as compared to the ones of an isolated quantum well. 

As a future perspective, the theory presented here can be generalized to explore the interesting and fascinating 
scenario of a time-modulated cavity. Even in the case of a vacuum input for the photonic and electronic polarization 
fields, a non-adiabatic temporal variation of the vacuum Rabi frequency (which is possible in the case of intersubband 
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transitions in doped quantum welli) is expected to produce an output radiation of correlated photons^, an effect rem- 
iniscent of the dynamical Casimir effect. The calculation of the quantum vacuum radiation spectra in the modulated 
case will be possible using and generalizing the comprehensive formalism developed in this paper. 
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APPENDIX A: GENERALIZATION OF THE INPUT-OUTPUT THEORY INCLUDING ALSO THE 
ANTI-RESONANT COUPLING TERMS IN THE BATH HAMILTONIAN 



In this Appendix, we generalize the input-output theory of the present system in order to include also the anti- 
resonant terms in the coupling with the photonic and electronic baths. These additional Hamiltonian terms make 
the calculations more cumbersome, but the results are still analytical. We have verified explicitly that when the 
broadening of the cavity and electronic transition modes is moderate, the effect of the anti-resonant terms in the bath 
Hamiltonians are negligible. The bath Haniiltonians including the antiresonant terms are 

k 

and 

k 

The generalized damping kernels 



and the generalized Langevin forces 



7i2,k(i) = ri2,k(t) - rl2^_kW > 

/cav,k(i) = Fcav,k{t) - ^cav,-k(*) ' 



(A3) 
(A4) 
(A5) 



The quantum Langevin equations in frequency space become 



(A6) 



where 



Mi:: 



2Dk -UJ- I7cav,k(w) 
-2Dk + j7cav,k(w) 



( ?k(w) ^^ 
6k(w) 
atk(-'^) 



( /cav,k(w) ^ 

-/cav,k(ti') 
V -/l2,k(c^) / 



, 



L012-UJ ~_ 1712, k(w) 
1712, k(w) 



2Dfc - i7cav,k('^) 

~i^R,k 

-Wcav.A; " SDfe ~~ UJ + «7cav,k(w) 
~i^R,k 



(A7) 



— 'i^R.k 

-j7i2,k(w) 

i^R,k 



-tJl2 - LJ + 1712, k(t^) / 

(A8) 
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the output operators q;°^ and f3°'^ are related to the input and cavity ones by: 



oout 



am 
Pq.k 



(AlO) 



As previously, Eq. I|A7I) can be solved exactly: 



/ ak(w) \ 

V &Lk(-^) / 



/ /cav,k(t^) \ 
/l2,k(w) 
-/cav.k('J-') 
V -/l2,k(t^) / 



(All) 



where 



k,cj J 



The relation between the input and output operators is 



(A12) 



g,k 



g',k 



(A13) 



with the generalized matrix 



UfriK u;) = l- 25R{f eav.k(^)} [GfriK CO) + Girik, u) - afr(k, CO) - QlPik, CO)] , (AM) 

ph -A- 

ufriKco) = -25R{fi2,k(c.)}^ [GfTi^^co) + GiriKco) - ar(k,c^) - gin^^^)] , (ais) 

N'.k 

z:?ir(k, c.) = _25R{feav,k(c.)}^ [^fr(k, ..) + ^"ir(k, - Girik, u^) ~ gin^, u)] , (ai6) 

<,k 

(k, c.) = 1 - 25R{fi2,k(c^)} [^"ir(k, ^) + ^4T(k, c^) - ^17 (k, uj) - ^f^ (k, c.)] . (A17) 

(A18) 

We have verified explicitly that the matrix W^'^"'(k, oj) is unitary. For the parameters used in this paper, iY^'^"'(k, lo) ~ 
i?(k, w), where Z7(k, w) is the unitary matrix (|38I41|) in the absence of the bath antiresonant terms. 
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Note that the anti-resonant terms Hanu in the Hamiltonian are due to the same vacuum Rabi coupling responsible for 
the resonant terms Hres. Hence, this problem, although there is some formal reminiscence, is different from the case of 
antiresonant transitions and Bloch-Siegert-like shifta— induced by a classical Rabi coupling produced by an intense and 
coherent pump field, which is comprehensively studied in Refl23. 



